Time-dependent transport via a quantum shuttle 
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We present a theoretical study of time-dependent transport via a quantum shuttle within the 
non-equilibrium Green's function technique. An arbitrary voltage is applied to the tunnel junction 
and electrons in the leads are considered to be at zero temperature. The transient and the steady 
state behavior of the system is considered here in order to explore the quantum dynamics of the 
shuttle device as a function of time and applied bias. The properties of the phonon distribution of 
the oscillating dot coupled to the electrons are investigated using a non-perturbative approach. We 
derive a relation for the oscillator momentum charge density correlation function which is an inter- 
esting physical example for the visualization of shuttling phenomenon. We consider the crossover 
between the tunneling and shuttling regimes for different values of the key parameters as a function 
of applied bias and time. We also consider the energy transferred from the electrons to the oscillat- 
ing dot as a function of time. This will provide useful insight for the design of experiments aimed 
at studying the quantum behavior of a shuttling device. 
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I. INTRODUCTION 

In recent years, there has been a great interest in 
the study of transport phenomena through nanoelec- 
tromechanical systems (NEMS)[l-4] both at an applied 
and at fundamental level. This attention has been fo- 
cused due to the fact that novel transport regimes arise 
due to strong interplay between electrical and mechan- 
ical degrees of freedom, which promise both the con- 
cept and realization of these systems as a new gener- 
ation of quantum electromechanical devices. A large 
number of new experimental techniques have been devel- 
oped to fabricate and perform experiments with NEMS 
in the quantum regime; Examples of high-frequency me- 
chanical nano-structures that have been produced are 
nano-scale resonators [ a |6 | , semiconductor quantum dots 
and single molecules [7Hll| . cantilevers [TH E3], vibrating 
crystal beams [1] and more recently graphene sheets 14 1 
and carbon nanotubes[l5l. fl(| . Anticipated quantum 
effects [1, H, [H, EH Il7| include observation of single 
phonons as well as ground state effects such as zero point 
displacement, ultra small mass sensing, measurement of 
extremely weak forces and suppression of scattering in 
coupled electromechanical systems. The strong interplay 
between a phonon mode and electrons is crucial for ob- 
serving the quantum effects in the transport of charge 
through electromechanical systems. 

A special kind of electromechanical system where the 
interplay between the electrical and mechanical degrees 
of freedom drastically ch ang es the transport properties is 
a shuttling device @, EH, EMI- The characteristic compo- 
nent that gives the name to these devices is an oscillating 
quantum dot of nanometer size that transfers electrons 
one-by-one between a source and a drain lead. A shut- 
tling device is a particular example of a NEMS device, 
an exciting generation of electronic devices such as: sen- 
sitive mechanical charge detector, logic gates, sensors, 
and memories. Beside these applications, a charge shut- 
tle may increase the understanding of quantum transport 
properties of NEMS and related quantum phenomena in 



nanoscopic systems. In an original proposal for a charge 
shuttle by Gorelik, et al[2(|., the investigation of the 
quantization of charge in the strong Coulomb blockade 
regime was discussed. This system consists of a movable, 
micro-size metallic grain spatially confined in a harmonic 
potential between source and drain leads. In the strong 
Coulomb blockade regime only a single excess electron 
at a time is allowed to occupy the grain. When a suf- 
ficiently high bias is applied between the leads, a sin- 
gle electron may tunnel onto the grain, and due to the 
electrostatic field between the leads, the charged grain 
is driven towards the drain, where the electron tunnels 
off. Due to the harmonic potential the uncharged grain 
is forced back towards the source and the process is re- 
peated. A mesoscopic shuttle system similar in some 
respects to that proposed by Gorelik et al. has been 
experimentally realized fl2l l2lj. This situation has been 
investigated both in the incoherent and in the quantum 
case. The shuttling phenomenon has been observed in 
Parks et al., experiment |7j, where a Cgo molecule is os- 
cillating between two leads. In this experiment, Parks et 
al., have observed the quantization of current for various 
applied voltages, which exhibits step like features within 
the current- voltage (I-V) curve. The molecular shuttling 
between the electrodes in this experimental setup has 
been done at low amplitude and high frequency but pro- 
vided a key evidence of the involvement of phonon levels 
in changing the current properties. This quantized be- 
havior in the conductance characteristics of NEMS de- 
vices has been observed in a similar experiment |22| . An- 
other method to drive the oscillations is to use an external 
AC electric field acting on a cantilever (Erbe et al..)[12j. 
where the amplitude can be tuned independently of the 
source-drain voltage bias. More recently [l9j. I-V charac- 
teristics of a quantum shuttle device have been observed 
at room temperature. These experiments are very close 
to the original proposal of Gorelik, et al. 

In general, there are two different theoretical formu- 
lations that can be used to study the quantum trans- 
port in nanoscopic systems under voltage bias. Firstly, 
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a generalized quantum master equation approach[23T- 
I2II and secondly, the nonequilibrium Green's function 
(NEGF)[13,[2l formulation. The former leads to a sim- 
ple rate equation, where the coupling between the dot 
and the leads is considered as a weak perturbation and 
the clcctron-phonon interaction is also considered very 
weak. In the latter case one can consider strong leads 
to system and electron-phonon coupling. The nonequi- 
librium Green's function technique is able to deal with a 
very broad variety of physical situations related to shut- 
tling transport 30] at molecular[3l| levels. It can deal 
with strong non-equilibrium situations and very small 
to very large applied bias. In the early seventies, the 
nonequilibrium Green's function approach was applied 
to mesoscopic j33 - l35j | transport by Caroli et al., where 
they were mainly interested in inelastic transport effects 
in tunneling through oxide barriers. This approach was 
formulated in an elegant way[28|, |36U39| by Mier et al, 
where they have shown an exact time dependent expres- 
sion for the non-equilibrium current through mesoscopic 
systems. In this model an interacting and non-interacting 
mesoscopic system was placed between two large semi- 
infinite leads. Furthermore, NEGF has been applied in 
the study of shot noise in chain models [4(| and disordered 
junctions[4l[ while noise in Coulomb blockade Joseph- 
son junctions has been discussed within a phase corre- 
lation theory approach 42] . In the case of an inelastic 
resonant tunneling structure, in which strong electron- 
phonon coupling is often considered, a very strong source- 
drain voltage is expected for which coherent electron 
transport in molecular devices has been considered by 
some workers [i^ within the scattering theory approach. 
Recently, phonon assisted resonant tunneling conduc- 
tance has been discussed within the NEGF technique 
at zero temperature[44|. To the best of our knowl- 
edge, in all these studies, time-dependent transport prop- 
erties via a quantum shuttle have not been discussed 
so far. The development of time-dependent quantum 
transport for the treatment of nonequilibrium system 
with phononic as well as Fermionic degrees of freedom 
has remained a challenge since the 1980's[13, US l36l — 
l39l HH . Generally, time-dependent transport properties 
of mesoscopic systems without an oscillating dot have 
been reported[46| and, in particular, sudden joining of 
the leads with a quantum dot or molecule have been 
investigated [45l |47| for the case of a noninter acting quan- 
tum dot and for a weakly Coulomb interacting molecu- 
lar system. Strongly interacting systems in the Kondo 
regime have been investigated (48j,|49(. More recently |50j. 
the transient effects occurring in a molecular quantum 
dot described by an Anderson-Holstein Hamiltonian has 
been discussed. 

In most of the theoretical work on shuttle based de- 
vices since the original proposal [2(|, mechanical degree of 
freedom has been described classically/semiclassically or 
quantum mechanically using the quantum master or rate 
equation approachj20l [23l - [26| . In the original proposal, 
the mechanical part was also treated classically, including 



the damped oscillator, and assuming an incoherent elec- 
tron tunneling process. This approach is based on a per- 
turbation, weak coupling and large applied bias approx- 
imations, whereas the Keldysh nonequilibrium Green's 
function formulation can treat the system-leads and 
electron-phonon coupling with strong interactions (5lj for 
both small and large applied bias voltage. Moreover, the 
theories of the Jauho group [24] and Armour 23] fail to ex- 
plain the low bias regime. The transport properties have 
been described and discussed semi-classically/classically 
but need a complete quantum mechanical description. 
A theory beyond these cases is required in order to fur- 
ther refine experiments on quantum shuttle based NEMS 
devices. In the quantum transport properties of a quan- 
tum shuttle; the quantized current can be determined 
by the frequency of the quantum mechanical oscillator, 
the interplay between the time scales of the electronic 
and mechanical degrees of freedom, and the suppression 
of stochastic tunneling events due to matching of the 
Fermionic and oscillator properties. For this reason, a 
first complete analytical and qualitative description of 
the quantum transport properties of a simple quantum 
shuttle device is performed and discussed in detail in this 
work. 

In the present work, we shall investigate the time evo- 
lution of an oscillating quantum dot as a reaction to a 
sudden joining to the leads. We employ the nonequi- 
librium Green's function method in order to discuss the 
transport properties of a quantum shuttle device. In this 
system, we include the time-dependent hopping between 
the oscillating dot and the leads to enable us to con- 
nect the leads to the dot at a finite time. An advantage 
of the time-dependent non-equilibrium Green's function 
approach considered in this work is that it eliminates the 
one major criticism of the Keldysh approach, which is 
the lack of any clear initial conditions. We assume strong 
dissipation of the oscillating dot, where we consider the 
average over single particle properties, implicitly assum- 
ing strong damping of the mechanical oscillator. Simply 
averaging over a single particle expression is a valid de- 
scription of a many particle problem, as long as any exci- 
tation caused by one electron has been dissipated before 
the next electron arrives. This is the strong dissipation 
regime. We describe the electronic state of the dot as 
a two-state system (empty/occupied). The electrons in 
the leads are considered to be at zero temperature. We 
consider an arbitrary finite chemical potential difference 
between the right and left leads. The transient and the 
steady state behavior of the system is considered here 
in order to explore the quantum dynamics of the shuttle 
device as a function of time. This is a fully quantum 
mechanical formulation whose basic approximations are 
very transparent, as the technique has already been used 
to study transport properties in a wide range of meso- 
scopic systems. In most of the existing literature very 
big chemical potential difference is considered while we 
include a range from very small to very large. In our cal- 
culation inclusion of the oscillator is not perturbative. As 
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the recent studies are beyond the perturbation the- 
ory, a non-perturbative approach is required beyond the 
quantum master equation or linear response. Hence, our 
work provides an exact analytical solution to the current- 
voltage, correlation functions, average energy, coupling of 
leads with the system, very small chemical potential dif- 
ference and includes both the right and left Fermi level 
response regimes. For simplicity, we used the wide-band 
approximation, where the density of states in the leads 
and hence the coupling between the leads and the dot was 
taken to be independent of energy. Although the method 
we used does not rely on this approximation. This pro- 
vides a way to perform transient transport calculations 
from first principles while retaining the essential physics 
of the electronic structure of the dot and the leads. An- 
other advantage of this method is that it treats the in- 
finitely extended reservoirs in an exact way in the present 
system, which may give a better understanding of the es- 
sential features of charge shuttle based NEMS device in 
a more appropriate quantum mechanical picture. 



II. FORMULATION 

The model consists of a moveable quantum dot sus- 
pended between the right and the left semi-infinite leads. 
This model is a close analogue to the dot attached to 
the tip of a cantilever or connected to the leads by some 
soft material or embedded into an elastic matrix. In this 
model the center of mass of the nanoparticle is confined 
to a harmonic potential. In the most recent experimental 
realizations of a quantum dot shuttle (QDS)[1_9( at room 
temperature, due to its small diameter, the quantum 
dot has a very small capacitance and thus has a charg- 
ing energy that exceeds the thermal energy k^T. Due 
to Coulomb blockade we consider only one excess elec- 
tron can occupy the device. We consider the electronic 
state of the central quantum dot as empty or charged 
(a two level system). The tunneling amplitude of the 
electrons depends exponentially on the position of the 
central island and thus electrons can tunnel between the 
leads and the quantum dot. This is due to the exponen- 
tially decreasing/increasing overlapping of the electronic 
wave functions. 

The coupled system is described by a single electronic 
level of energy eo and a nanomechanical oscillator with 
frequency oj and mass \x. An electrostatic force (eE) acts 
on the mechanical grain when the quantum dot is charged 
and gives the electrical influence on the nanomechanical 
dynamics of the system. The electric field E is generated 
by the voltage drop between the left and the right leads. 
In our model, though, it is kept as an external parameter. 
The quantum dynamics of the nanomechanical oscillator 
is represented in terms of position and momentum op- 
erators. In terms of creation and annihilation operators 
for the nanomechanical oscillator excitations the Hamil- 
tonian of our simple system is@, HH-IH, H El El 



Hq — Hdot-ph + H[ e 
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#dot-ph = [eo + r)(tf + b)] ejeo + H ph . 
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= fiw(&t&+ 1), 



(1) 



(2) 



(3) 



where eo is the single energy level of electrons on the dot 
with c ,Cq the corresponding creation and annihilation 
operators, and an oscillator of frequency u>, mass /j, and 
b' , b are the raising and lowering operator of the phonons. 
The parameter rj = physically represents an effective 
electric field in the capacitor seen by the moving dot be- 
tween the leads, which we shall refer to as the coupling 
strength between the moving dot and the electrons on 
the dot given as = eEl, where e is the charge of elec- 
tron, E is the strength of electric field and I = A /— is 

the zero point amplitude of the oscillator. The remaining 
elements of the Hamiltonian are 



H 



leads 



fc 



£ j C j C j' 



(4) 



H 



leads-dot 



= AH a = 



N 



^y a (t)(c* C „ + 4), (5) 



where N is the total number of states in the lead, a = 
L,R, and j represents the channels in one of the leads. 
We include the time-dependent hopping V a (t) to enable 
us to connect the leads a to the moving dot at a finite 
time. An advantage of this approach (joining of leads 
to the dot at t=0) is that it eliminates the one major 
criticism of Keldysh approach, namely the lack of any 
clear initial conditions. V a (t) is written as 

V a (t)=V (t)e*& (6) 
where £ — ^ is the inverse tunneling length (£), — 
stands for Vt(i) = Vo{t)e~^ x and + stands for Vn(t) — 

Vo(t)e + ^ x . For the time-dependent dynamics, we shall 
focus on sudden joining of the leads to the moving dot 
at t = 0, which means Vo(t) = V6(t), where 9{t) is the 
Heaviside unit step function. The displacement operator 
is given as 



a i(6t + b) 

x = F — 

V2 



(7) 



The total Hamiltonian of the system is thus H = Hq + 
AH a . We write the eigenvalues and the eigenfunctions 
of Hdot-ph as 



e + tuj(n + 5 ) 



(8) 
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^ m {K,XQ ^ 0) = A ro exp[-^]H m (^)exp[-ii^a;o] 

(9) 



^ n (K,x = 0) = A n exp[- 



1 



(10) 



for the occupied, xq ^ and unoccupied, 



respectively, where A n 



dot 



2^ 2 



.I'O 



A 



is the displaced oscillator equi- 
librium position due to the coupling to the electrons on 
the dot, and H n (lK) are the usual Hermite polynomials. 
Here we have used the fact that the harmonic oscilla- 
tor eigenfunctions have the same form in both real and 
Fourier space. 

In order to transform between the representations for 
the occupied and unoccupied dot we require the ma- 
trix with elements $„ !?n = / \l>* (K, x = 0)^ m (K,x ^ 
0) dK, which may be simplified [54] as 



'2l™-"l minfn,ml! , -, 9X 
—exp(-±x z ) 



max[n, m]\ 



\m—n\ j^\m-n\ 



(Kh (11) 



where x = ^f- is the ratio of equilibrium position of the 
displaced oscillator to the zero point amplitude of the os- 
cillator and L™~ n (x) are the associated Laguerre polyno- 
mials. The position of the resonant level with respect to 
the chemical potential in the leads is thus affected by x of 
the nanomechanical oscillator, which in turn affects the 
transport properties of the junction through the device. 



III. THE SELF-ENERGY AND THE GREEN'S 
FUNCTION 



In order to calculate the analytical solutions and to 
discuss the numerical results of the transient and steady 
state quantum dynamics of the nanomechanical shuttle 
device, our focus in this section is to derive an analyti- 
cal relation for the time-dependent effective self energy 
and the corresponding Green's function. The effective 
self-energy represents the contribution to the moving dot 
energy, due to interactions between the oscillating dot 
and the leads it is coupled to. In obtaining these results 
we use the wide-band approximation only for the sim- 
plicity, although the method we are using does not rely 
on this approximation. The retarded self-energy of the 
oscillating dot due to each lead is given bv(28l. I36U391 1451] 
(see Appendix) 



S r (Mi 



ir. 



^ 2 )^-^)S r ( T ) (12) 



with the matrix £ r (T-) 

= exp[(7T x) 
x (T7) K " 




2l™o-™o| 



min[rio, n' }\ 



mm[no ,n 



max[no, n ] 



(13) 



where £ a (i,ii) = (S r (i,ti))* with a representing the L 
or R leads, the dimensionless tunneling length 7 = 



^2-, and the 



signs stands for the left and the 



right leads respectively. In the present representation, 
the effective self energy is given in terms of a matrix 
due to the position dependence of the tunneling matrix 
elements. 

The final result for the lesser self energy matrix may 
be written as 



V < (t 1 ,t 2 )=iT a 6(t 1 )e(t 2 ) x £<(t) 



(14) 



-™ 27T 



exp[-k Q (t! - t 2 )], 



where the matrix element of the lesser self-energy is writ- 
ten as 



V<(T)[V<(T)] T , 



(15) 



cxp[- 



[T(4 £ )] n 



S<(T) 

where V<(t) = V<(t) 

This model represent the interplay between two physi- 
cal time scales of the system, i.e., the oscillator frequency 
(ui), and the tunneling rate (T a ). The model also shows 
very interesting interplay between three physical length 
scales of the system, i.e., dimensionless tunneling length 
(7), zero point amplitude (I) and the displaced oscillator 
equilibrium position to the zero point amplitude of the 
oscillator (x), which are actually affected by the weak and 
strong coupling dynamics and small and large tunneling 
rate. 

We solve Dyson's equation using i?dot-icadsj as a pertur- 
bation in terms of matrix representation. In the presence 
of the oscillator, the matrix of the retarded, G r (t, t\) and 
advanced, G a (t2,t') Green's functions on the dot, with 
the phonon states may be written as 

G r (t, h) = -\0{t - ti) cxp[-i(M)(t - ti)], h>0 

(16) 

where G r and M are matrices in nanomechanical oscilla- 
tor space with indices m, no, and the matrix M may be 
written as 

M = (E m S m . na + no ), E m = e Q + (m + |)fiw - A, 

(17) 

G a (t 2 ,t r ) = (G r (t / ,t 2 ))* 

= +id(t' - t 2 ) exp[-i(M*)(t 2 - t% t 2 >0 
= G a n , k {t 2 ,t') (18) 
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The above Eqs. (12), (14), (16), and (18) will be the 
starting point of our examination of the time-dependent 
response of the coupled system. These functions are the 
essential ingredients for theoretical consideration of such 
diverse problems as low and high voltage, coupling of 
electron and phonons, transient and steady state phe- 
nomena. 



IV. THE DENSITY MATRIX p n<n , (t, t) AND THE 
CORRELATION FUNCTION 



The density matrix is related to the lesser Green's func- 
tion through Pn,n'{t,t) — —iG^ n ,(t,t') at t' = t, where 
the G5 n , (*,*') is the lesser Green's function [13, [H, 0- 
Ea| on the oscillating dot including all the contribu- 
tion from the leads. When t and t' > the lesser Green's 
function for the dot in the presence of the nanomechan- 
ical oscillator on the dot, with phonon states in the rep- 
resentation of the unoccupied dot, may be written as 



In terms of eigenvalues and eigenfunctions, the advanced 
Green's function is written as 



(24) 



with eigenvectors U* and its inverse such that U*TJt = I, 
and 



g a (t 2 ,t') = g?(t 2 ,t') 

= +i9(t' - t 2 ) exp[-i(e. i + iQ)(h - t')]8 id , 

(25) 

The lesser self-energy, 'E < (ti,t 2 ), in matrix space of the 
oscillator was derived in the previous section. 

The lesser Green's function can be calculated by using 
Eqs. (14,21, & 24) in equation (19) at t = t' as 

= f f dt 1 di 2 #Ug r (t,i 1 )UTS < (t 1 ,i 2 )U*g a (t 2 ,i)U t $ t , 
Jo Jo 

= [ [ dM*2 y)[su] Tl , fl g r (t,f 1 )UT 

Jo Jo „ h „ 



g,b,ct 

xS<( T )S<(t 1 ,t 2 )U*g a (t 2 ,t)[*U] t M „ 
= ^[*U]„, 3 [UTS<(T)U*] 3M « b (t)[*U] t M „ (26) 



9,f>,a 

where the function X g b(t) ^ s written as 



G<„<(M') 

= J2 I I dhdh 

1 , Jo Jo 

= f f #G r (t,t 1 )S<(t 1 ,t 2 )G a (t 2 ,t')* t dtidi 2 (19) 
Jo Jo 

whereas for t and t' < 0, the G< n ,(t,t') is equal to zero. 
G n,n> & includes all the information for the nanome- whic h can be written as 
chanical oscillator and electronic leads of the system, and 
m,k,no,ni,n,n' are the oscillator indices. The lesser 
self-energy, aTll (ti,t 2 ), contains electronic and oscilla- 
tor contributions. The electronic contributions are non- 
zero only when t\ and t 2 > 0. The retarded Green's 
function, G r m nQ (t, ti), is represented as 



Xg <b (t)= Ahdt^HtM^hMKit^t) (27) 
Jo Jo 



Xg, b (t) = / / dhdt 2 exp[-i(e s - iQ(t - h)} 
10 Jo 




e Pa +±Hw d£a 
-™ 27T 



exp[-ie Q (<i - t 2 ) 



G r (t, t x ) = -i6(t - h) exp[-i(M)(t - t x )), (20) 

Using the eigenvalues and eigenfunctions of M, the re- 
tarded Green's function may be written as 



x exp[-i(e b + iCb)(<2 -*)], 



(28) 



G r (t,ti)=Ug r (t,ti)U T , 



(21) 



Although G r ( a \t, t') is non-zero for t < 0, it is never re- 
quired due to the way it combines with Sa (i? t'). 
By carrying out the time integrations, the resulting ex- 
pression is written as 



with eigenvectors U and its inverse such that UU T = I 
(not UU f = I), and 

g r {tM) = 9 r S,ti) 

= -i0(t-ti)exp[-i(ei-iCi)(t-ti)]*ij, (22) 

where £i,Ci are the real and imaginary parts of the ith 
eigenvalue of the matrix M. The advanced Green's func- 
tion, k (t 2l t), is written as 

G°(t 2 , t') = +i6(t' - t 2 ) exp[-i(M*)(t 2 - t% (23) 



ir„ 



1 



2tT 7-oo " ( £ Q - £fc ~ Kb)(£a - £g + Kg) 

X jl + exp[i(e 6 - eg + i(Cb + ( g ))t) 
- exp[-i(e Q - e b - i( b )t] 

-exp[i(e a ^e g + iC g )t}\ (29) 
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For all values of b, <?,and a, Xg b(t = 0) must be equal to 
zero. The integral over the energy in the above equation 
is carried out[5a|- The final result for the above expres- 
sion is written as 



1 



ilV 



2lT £ b - Eg + i(Cb + (g) 



}[Y g a b + Z«bl (30) 



where we have added the contribution from the right and 
the left leads, which can be written in terms of a as 



XfOt 



(l + exp[i(e b -e s +i(C b + C s ))i]) 
1 ln[(e Fce - e b ) 2 + C b 2 ] 



e g ) 2 + Q] 



2 ln[(e Fct - c g 

1 / £Fa ~ £b 



tan" 



tan 



C 9 



(31) 



and 

ryOi 



= exp[i(e b - e g + i(Cb + Cs))*] 
x (- Ei[it(e Fa -e b - Kb)} + Ei[-it(e Fa -£ g + Kg)}) 
+ Ei[i<(e FQ - e g + Kg)} - Ei[-i<(e Fa - s b - Kb)}, 

(32) 

with e Fa being the right and the left Fermi levels and 
Ei(a;) the exponential integral function. Special care is 
required in evaluating the Ei(a;) to choose the correct Rie- 
mann sheets in order to make sure that these functions 
are consistent with the initial conditions & (t = 0) = 
and are continuous functions of time and chemical po- 
tential for each value of b,g, and a. 

Now using equation (26), the dot population may be 
written as 

Pit) = 5>n,»(*> *) = E - iG «,«(*> *)' ( 33 ) 
n n 

Next we derive a relation for the correlation function 
which is an interesting physical quantity to see the shut- 
tling dynamics of the system. With the help of equation 

(26) the correlation function Tr< pG < (t,t) >, where 
p = ^p(&t — 5) i s the momentum operator of the os- 
cillator, may be written as 



Tr 



UflLO 



(6 f - b)G<(t,t) 



< 

n,n- 



i(*,*)}>. (34) 



which we expect to be finite in case of shuttling as the 
dot should be occupied for positive momentum while zero 
for unoccupied dot. As we note that the momentum cor- 
relation function depends on off-diagonal elements of the 
lesser Green's function which can only be finite when 
nanomechanical oscillator is in a mixed quantum state 
rather than in a simple eigenstate. 



V. TIME-DEPENDENT CURRENT FROM 
LEAD a 

The particle current I a into the interacting region from 
the lead is related to the expectation value of the time 
derivative of the number operator N a = Yl a j c aj c <xj> 



T.,=-e{±*) = %([x t B\) 



and the final result for the current through each of the 
leads is written as (See Ref. 51 for detail) as 



I a (t) 



dti 

x Tr{G r (t,ti)H < (t 1 ,t) + G < (t,t 1 )-E a (t 1 ,t) 
- S r (i,t 1 )G<(< 1 ,i) - S < (t,i 1 )G a (t 1 ,t)}, 

j:(Z? + Z?) (35) 



where Z" is defined as 

Z?= f dtiTr{G r (i,ii)E < (t 1 ,0-S < (t,t 1 )G a (i 1 ,i)} 
Jo 

= [ d< 1 Tr{Ug r (t,t 1 )UTS<( T )S<(ti,t) 
Jo 

-S<( T )S<(t,t 1 )U*g a (t 1 ,0Ut} 

= ]T{[UTS<( T )U] 9 , sX ^) 

g 



(36) 



where Xq(t) is defined as 



Jo 



/ d*iexp[-i(e fl -iC fl )(t-ti)] 
Jo 



-i0(t-tx) 
' 2^ 

ds a exp[— iE a {ti — t) 

-co 

ir„ r FaHhu d£ 1 - exp[-i(g a - e g + Kg)t] 

27T ' £ a - € g + Kg 

iT 



2tt 



j{ln[(e Fct - Eg) + Kg} ~ Ei[it(e Fa - £ g + Kg)}} 

(37) 



Hence, the final result for Z" is written as 
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Z i a = S{E[ UTS< (T)Uk 9 

x { i ln[(e Fa - e g ) + iC 9 ] - Ei[i*(e Fa - e g + i( g )} } 

(38) 



and Z£ is written as 



Z£ = J dtiTr{G<(t,ti)S a (ti,t) - S r (t,t 1 )G<(ti,t)} 
= J d< 1 Tr|G < (t,t 1 )S a ( T )^^(ii)<5(i - t x ) 
+ E r (T)^e(ti)*(t-ti)G<(*i,t) 



Tr<j ^G<(i,t)S a ( T ) + ^E r ( T )G<(f,t) 



i£o 

2 - \-,-y- v . / ■ 2 

= ir Q / / d; 2 dt 3 

Jo Jo 

Tr{C"(t, t 3 )S<(T)S< (t 3) i 2 )G a (i 2 , t)£ a ( T )} 

= iT Q / / di 2 dt 3 
Jo Jo 

Tr{Ug r (i,t 3 )UTS<(T)S<(i 3 ,i 2 )U*g a (t 2 ,t)UtS a (T)} 
= ir Q {^[[UtS Q (T)U] M [UTS<( T )U*] 3M « b (t)}, 

5, 6, a 

(39) 

where Xg (,(*) is as in ([3"0]l. 

Now using Eqs. (001 E3 & [3S]) in equation (05]), we 
arrive at the final result for the current through lead a 

as 



/ Q (i) = -(Z« + Z 2 L + Z«), 



(40) 



where in calculating the left current we need and both 
the contributions Z^ and Z^ and for the right current Z^ 
is replaced by Zf-. As before, special care is required in 
evaluating the ln(x) and Ei(x) to choose the correct Rie- 
mann sheets in order to make sure that these functions 
are consistent with the initial conditions I a (t) = and 
are continuous functions of time and chemical potential. 



VI. AVERAGE ENERGY 

To calculate the energy transferred from the electrons 
to the nanomcchanical oscillator, we return to the density 
matrix p n ,n(t, t) given in Eq. (33). We may use the lesser 
Green's function or density matrix to calculate the energy 
transferred to the oscillator as 



Eph =< n > fko = 



(41) 



where the average evaluated using the diagonal element 
of the density matrix on the oscillating quantum dot. 
Note that the normalization in equation (42) is required 
as the bare density matrix contains both electronic and 
oscillator contributions. The trace eliminates the oscilla- 
tor part, leaving the electronic part. 



VII. DISCUSSION OF RESULTS 

The net (7t(t) — In(t)) average current through the 

system, the correlation function < pG < (t,t) >, total 
(ii(t) + In(t)) average current into the system, and the 
average energy of an oscillating dot between the leads 
are shown graphically as a function of time (2ir/u)) for 
different values of tunneling length, tunneling rate, and 
voltage bias. The following parameters were employed: 
the single energy level of the dot e = 0.5 and ep_R = 0. 
These parameters will remain fixed for all further discus- 
sions and have same dimension as fku. We are interested 
in small and large values of tunneling rates (T L / R ) from 
the leads, different values of the tunneling length (£) be- 
tween the oscillating dot and the electrons, and of the 
left chemical potential < e-pL < 2. The oscillating dot 
induced resonance effects are clearly visible in the numer- 
ical results. The tunneling rate of the electrons between 
the leads and the dot is considered to be asymmetric 
{Tr = r^expf— jx]) and we assume that the leads have 
constant density of states. The values of the tunneling 
rates has been chosen such that when the dot is occu- 
pied and the zero point has moved to x , these values are 
reversed. 



1.12ft 



1.10 



1.08 - 



1.06 



1.04 h 

j i i i i i i i i i_ 

20 

Time (In/at) 

FIG. 1: Time-dependent net average current 27r(ii(t) — 
In(t))/efko against time (2-k/ui) for fixed values of the param- 
eters: e = 0.5, €fr = 0, efl = 2, x = 0.21, fku — 0.63, Fl = 
foj/27r, r_R = r_L exp[— yx]) and 7 = 0.5. The electron per cy- 
cle behavior in this figure corresponds to the shuttling mech- 
anism. Units: all the parameters have same dimension as of 
fius. 

In Fig. Q] & [5] we have shown the net average current 
1-k{Il(S) — Iu(t))/ehuj flowing through the system as a 
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2n (7 z (t)-7 8 (t))/e* ui 




Applied voltage 



FIG. 2: Net average current (2ir(I L (t) — I R (t))/ehcj) flowing 
through the system as a function of both time (2n/u}) and of 
the left Fermi level for same values of the parameters as in 

fig- ID 



0.105 - 

3 

-? 0.100 - 



0.095 - 




10 20 30 40 50 

Time (2jt/co) 

FIG. 3: Net average current (2n(lL(t) — In(t))/eHu)) flowing 
through the system as a function of time (2n/us) for fixed 
values of the parameters in case of the tunneling mecha- 
nism: eo = 0.5, efr = 0, cfl = 2, Yiw = 0.63, Tl = = 
0.17iu;/27r, xo = 0.1? and 7 = 0.001. Units: all the parameters 
have same dimension as of Hu>. 

function of time (2n/oj) for fixed values of the Fermi 
level, £fl = 2, e^R = 0, tunneling energy, Tl — huj/2ir, 
the characteristic energy of the oscillator hu = 0.63, 
Tr = Tl exp[— 72;], xq — 0.3Z and finite 7 = 0.5. 
This (7 = 0.5) is the physical condition for shuttling 
in the quantum regime: when the zero point amplitude 
of the nanomechanical oscillator is half of the tunnel- 
ing length (C = In this figure, in which the os- 
cillator period is comparable with the left bare tunnel- 
ing rate (Tl = fiw/27r), the average current saturates at 
one electron per mechanical cycle (corresponding to av- 
erage current 2n(lL{t) — In(t)) / ehaj = 1) as the electrons 
are shuttled one by one from the source to the drain by 
the oscillating dot. We get the following physical pic- 
ture: every time an electron jumps onto the dot when 
= hui/2ir, the dot is near the source lead and subject 
to the electrostatic force eE that accelerates it towards 
the drain lead. Energy is pumped into the nanomechan- 
ical oscillator and the dot starts to oscillate between the 
leads. When the right tunneling time is high compared to 
the left tunnelling time, the oscillator dissipates this en- 
ergy into the environment before the next event occurs. 



This continuously drives the oscillator away from equi- 
librium and a stationary state is reached only when the 
energy pumped per cycle into the system is dissipated 
during the same cycle in the environment. With this 
setup, the system will be in the shuttle regime. We note 
that the frequency of the short period weak oscillations 
is (| cfl — eo|) and these oscillations are present even in 
the limit of small but finite 7. We conclude that this is 
a purely electronic process (plasmon oscillations). It is 
clear from the figure that in the strong tunneling case, it 
contains two beating frequencies, therefore we interpret 
this as due to a mixture of electronic and mechanical 
frequencies. In Fig. [2] the long time behavior of the net 
average current has steps as function of the left Fermi en- 
ergy which become more pronounced with increasing 7. 
In contrast, when both the bare tunneling rates are equal 
and smaller than the mechanical frequency, and 7 << 1, 
the dynamics of the system is similar to a double barrier 
resonant tunneling system, since the dot is static and 
far from both the leads. The tunneling dynamics of the 
net average current is shown in Fig. [3] for fixed values of 
Tl = = 0.lhuj/2ir, xo = 0.11 and the oscillator energy 
hu> =0.63. This confirms the tunneling behavior of the 
system oscillating with beating frequency which is due 
to mixture of electronic and nanomechanical degrees of 
freedom and is consistent with previous studies [5 1[. 
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0.00 
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FIG. 4: The momentum correlation function as a function 
of time (2tt/uj) showing the shuttling (solid line) and the 
tunneling (dotted line) behavior. Parameters for the shut- 
tling are eo = 0.5, fiu = 0.63, Tl = fkd/2n,xo = 0.21, 
Fr — r_Lexp[— yx]) and 7 = 0.5, and for the tunneling 
are eo = 0.5, €fr = 0, £fl = 2, Hu) = 0.63, Tl = Tr = 
O.lfkj/2-K, xo = 0. li and 7 = 0.001. All the parameters have 
same dimension as Hu>. 

The correlation function is shown in Fig. [4] & [5l as a 
function of time {2tt /ui) and left Fermi energy for fixed 
values of tunneling length £ = 21, xq — 0.51 , frw = 0.63 
and tunneling rates Tl = fojj/2ir, Tr = Tl exp[— jx]. 
We distinguish between the shuttling and tunneling dy- 
namics of the system for long time for different tunnel- 
ing rates and 7: For shuttling (solid line) Tl = fooj/2ir, 
Tr = Tl exp[— 7a;], 7 = 0.5, and for tunneling (dotted 
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FIG. 5: The momentum correlation function as a function 
of time (27r/w) and left Fermi energy showing the shuttling 
behavior. Parameters for the shuttling are eo = 0.5, fhj = 
0.63, T L = fnu/2-K,x = 0.21, F R = r L exp[- 7 a;]) and 7 = 1. 
All the parameters have same dimension as of tku. 



line) T L = T R = 0.1fiw/27r, x = 0.1/ , and 7 = 0.001. 
Why have we found this particular type of structure? 
We know that, in the ideal shuttling regime, we ex- 
pect the dot to be occupied when the momentum is 
positive and unoccupied when the momentum is neg- 
ative. Hence, the correlation function is positive. In 
contrast, when the transport is independent of the me- 
chanical oscillator, the correlation function is zero as 

< pG<(t,t) >=< p >< G<(t,t) > = 0. We observe 
that the long time value of the correlation function has 
steps as a function of the left Fermi level as shown in 

fig- El 




50 100 

Time (2^/6j) 



200 



FIG. 6: Total average current (li,(t) + Iii(t)) flowing onto 
the dot as a function of time (2tt/lo) for fixed values of eo = 

0.5,6 FR = 0,€ FL = 2,Hu = 0.63, T L = F R = 0.lhuj/2TT,X = 

0.21 and 7 = 0.1. All the parameters have same dimension 
as tiu). This current (solid line) is equivalent to the rate of 
change of dot population ^p{t) (dashed line) as a function of 
time for same parameters as of current. In this figure, solid 
and dashed lines have same values at all points. Units: all 
the parameters have same dimension as of fej. 

Next, we have shown the total average current (Ii,(t) + 
In(t)) flowing onto the oscillating dot in Fig. H as a 
function of time {2tt/uj) for fixed values of Tl — Tr , — 
0.1fio;/27r, fku = 0.63, epi? = 0, £fl = 2, and finite value 
of the tunneling length 7 = 0.1. This current (solid line) 
is equivalent to the rate of change of the dot population 



(dashed line) for the same parameters. In this figure, 
we can not distinguish the solid and the dashed line. 
This confirms that our analytical results are consistent 
with the equation of continuity, Jz(i) + In,(t) = j|p(t), 
and hence, with the conservation laws for all parameters. 
We interpret the long period oscillations as a function of 
time (2ir/u)) is the mechanical frequency of the system 
and short period weak oscillations with the electronic 
frequency (|£fl — eo]), which is consistent with the net 
average current results as well. 




FIG. 7: Average energy transferred to the oscillator as a 
function of time (2tc/lo) for fixed values of ttoj = 0.63, Fl = 
tko/2n,xo = 0.5Z, F_r = r_Lexp[— 72;]) and 7=1. Units: all 
the parameters have same dimension as of fej. 




Applied voltage 



Time (2n/u>) 



FIG. 8: Average energy transferred to the oscillator as a func- 
tion of time (2n/oj) and left Fermi energy for all the parame- 
ters same as in fig. [7] 

Finally we have shown the average energy of the 
nanomechanical oscillator as a function of time (2tt/uj) 
and of the left Fermi energy in Fig. [7] & [8] for fixed values 
of tunneling rates Tl — 0.1, Tr — Tl exp[— jx], cfr. = 0, 
xq = 0.31 and tunneling length 7 = 0.5. These param- 
eters corresponds to the shuttling regime of the system. 
We found constant average energy of the nanomechan- 
ical oscillator for long time. This constant average en- 
ergy increases with increasing Fermi level as shown in 
fig. [HI We consider two possible interpretations of this 
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structure. One is presented in the previous work: that 
the oscillator is in any of its pure states and therefore 
the potential seen by the electrons is time independent. 
Secondly, this is an average over several time-dependent 
processes with different phases such that any Rabi like 
features cancel. Hence, we note that only the second in- 
terpretation is consistent with the correlation function. 
For the case of the tunneling regime of the system, the 
average energy of the nanomechanical oscillator is shown 
in Fig. IH1 as a function of time (2tt/uj) for fixed values 
= Fr = 0.1Hui/2ir, fko = 0.63, ep_R = 0, ep^ = 2, and 
for very small value of the tunneling length I = 0.001C- 
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FIG. 9: Average energy transferred to the oscillator as a func- 
tion of time (2n/u>) for fixed values of eo = 0.5, €fr = 0, efl = 
1, fko = 0.63, F L = F R = 0.1fiw/27r, x = 0.51 and 7 = 0.001. 
Units: all the parameters have same dimension as of fko. 

In Fig. [71 [S] & HI h appears that the oscillator is not in 
its ground state at t=0. However, if we look at the struc- 
ture of the lesser self-energy, the quantity in the mid- 
dle refers to the state of the oscillator when the dot is 
empty and the oscillator is in its ground state. Whereas 
the whole self-energy refers to the state of the oscillator 
when an electron has passed through the barrier onto the 
dot. At this point the oscillator is not necessarily in its 
ground state. When the electron is in the leads, the os- 
cillator is in its ground state but when the electron is 
on the dot two things happen: the representation of the 
lesser self-energy is changed and the other is the effect 
of the exp[=F7a;] operator which kicks the oscillator out 
of its ground state. This comes back into a saturated 
state for long time with finite value of average energy. 
This could be a mixed or a thermal state. When taken 
with the finite correlation function for the same param- 
eters, we conclude that the long time average energy of 
the nanomechanical oscillator corresponds to the mixed 
state. 



VIII. SUMMARY 

In this work, we analyzed the time-dependent trans- 
port via a quantum shuttle by using the nonequilibrium 



Green's function approach without treating the electron 
phonon coupling as a perturbation. We have derived an 
expression for the full density matrix and discuss it in 
detail for different values of the tunneling length and the 
tunneling rate. Using the full density matrix calculation 
we have shown the oscillator momentum charge density 
correlation function which distinguishes the shuttling and 
the tunneling regime of the system. We derive an expres- 
sion for the current to see the effects of the coupling of the 
electrons to the oscillator on the dot and the tunneling 
rate of electrons to resolve the dynamics of the nanome- 
chanical oscillator. We found a positive correlation func- 
tion for the nanomechanical oscillator in the shuttling 
regime and zero in the case of the tunneling regime. The 
shuttling regime will occur when the electron left jump 
rate introduces continued kicks on the island. This hap- 
pens when the left jump rate is close to the oscillation 
frequency of the island — ftw/2ir. By setting an ap- 
propriate ejection rate (r^ = Tl cxp[— 72;]) and I = 0.5(, 
there exists a condition where the island will keep os- 
cillating between the leads. In contrast, in the limit of 
equal bare tunneling rates, smaller than the mechanical 
frequency, and when 7 < < 1 , the dynamics of the system 
is similar to a double barrier resonant tunneling system, 
where we found periodic oscillations as a function of time 
(27t/cj), which corresponds to the mechanical frequency 
of the system. Furthermore, we discuss the average en- 
ergy transferred to oscillator as a function of time. We 
have found a mixed or a thermal state shown in the aver- 
age energy of the nanomechanical oscillator for long time 
which is consistent with the finite correlation function for 
the same parameters. These results suggest further ex- 
periments for NEMS to go beyond the classical dynamics 
of a shuttle device. 



Appendix A: 

In this appendix we derive an analytical relation for 
the time-dependent self-energy The effective self-energy 
represents the contribution to the moving dot energy, due 
to interactions between the oscillating dot and the leads it 
is coupled to. The retarded self-energy of the oscillating 
dot due to lead a is given by 

2^„,(Ml) = 4>n,rn(t)g^ a {tMWn', m (tl), (Al) 

where the uncoupled Green's function g r aa {t,t\) in the 
leads is 



+00 



5a,a(Ml) 



3 

+00 



5a,,' (Ml) 



n a de a g r a {t,tx) (A2) 



i0(t — ii) / de a n a exp[—i(e a + (m + ^)fuv)(t — ti)], 
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and £ / Nn.de., where j stands for every ehan- 

J — oo 

nel in each lead, n a is the density of states in lead a, 
9a,j(^>*i) ^ s ^ ne uncoupled Green's function for every 
channel in the leads. The tunneling matrix <j) n ,m{i) is 
written as 



where x is the coordinate of the oscillator, £ being the 
inverse tunneling length. 



<t>n,m(t) = ( V (i)e : 



Ale ^^Hr, 



V (t)e 



A m e 2,2 Hr. 



X — Xq 
I 



(t) 



da;, 



Using the expressions of the tunneling matrix and the 
(A3) uncoupled retarded Green's function in equation (Al), 
the retarded self-energy may be written as 



+°° ,2 



K,n>(t,ti)=J2- i6 ( t - t ^ I Nds " y ^ n a exp[-i(s a + (m + $)hu){t - h)] 
m 

x / dxA n e — «a— H„ ( — j — J F e« x ^e li*" H m (yj 

x / d.T'A^e" 1 "^ H„/ (^^) foe** A^e"^ H m (y) , (A4) 



+ 00 

-iO(t - h) \V\ 2 0(i)0(*i)n Q y exp[-i(e Q + (m + - h)]de c 



^)n n (^)^exp(-f)n : (^ 

x | dx'A* n , exp (- (X, ~/° )2 ) Hn' V e^'A m exp (-^) H m (y) , (A5) 

= -in Q |F| 2 0(t)6(ti)9(t - h)2iv6(t - h) 
x £ / d^„ exp (-^=^) H„ (^) exp Hm (f ) 

x | da^, exp (- ^"r ^ ) H„> (^y^) V e^'A m exp (-^) H m (y) , (A6) 
where r a =4ir \ V\ 2 n a . Using the completeness identity, 

J2 A *n ex P (-^Jt) Hm (y) A m exp (~yy-) H m f y) = 5(x - x') 



the above equation can be simplified as 

Ky (t,t 1 ) = -^9(t 1 )5(t-t 1 ) (A7) 

x / d x'A n exp (~^-/^) H„ (^) exp (-^^) H n ' (^) , 

I 



After integrating the above expression, we arrive at the following final result for the retarded self-energy[5 

£L„'(Mi) = -^(iiW-MSm (A8) 
E;(Mi)E r (T) = S r (Mi 
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with 



(A9) 



where t\) — ti))* with a represents the L or 

R leads, and the matrices ip n ,n' is given as 



/2l n '-"l min[ra,n']! 



R ^(T7)'"-"'L!;; n[ ; i „, ] (-2y), 

(A10) 



where the dimensionless tunneling length 7 = 

and the — , + signs stands for the left and the right 
leads respectively. 

Similarly, the lesser self energy may be calculated as 

£<„,(*!, t 2 ) = ^,m(*l)<a(*l,*2)C', ro (*2), (All) 

where a (ti,t 2 ) is the uncoupled Green's function in 
the leads and is given as 

+00 

ffa,a(*i>*2) = ^S 5 ^* 1 '* 2 ) = y n a de a g<{t u t 2 ) 

3 — OO 

(A12) 



;/; = 



+00 

/ 



2in a /(e Q )-B m <fe a exp[-i£ Q (ti - i 2 )], 



We are going to consider that temperature is zero in the 
leads. So, we assume, when the oscillator appear in the 
leads it will always be in its ground state such that m = 
(B m = S m fi) and the Fermi distribution function can be 
replaced with the theta function which in turn limits the 
integral. With the help of these approximations, and 
using the expressions for the tunneling matrix and the 
uncoupled retarded Green's function into equation (All), 
the lesser self-energy may be written as 



/• £ f» + 5 ftw \V\ 2 

^,n'(*i>*2)= / N2in a de a l -^9(t 1 )e(t 2 )cxp[-ie a {t 1 -t 2 )} 

x / dxAn exp (-^#) H„ (^) Voc^A* exp (-£ 
x Jdx'A* n/ 

efe a exp[-ie a (ti-t 2 )] 

-OO 

x / i-,A, exp H„ ( '-^A) V^M exp < (X) 



exp H,„ % e«^oexp (-<g! ) , (A13) 



2/ 2 j ""■ v j / """^ v 2Z 2 



X |d^ 

= il^O^) / ^ exp[-i£ a (t! - t 2 )\ 

J-00 2lT 

x / d,A„ exp (-^f)!) H„ (^) V e^A* exp 

x /d^exp (-^^) H n ' (^) Foe^Aoexp (-^ ) , (A15) 
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where r Q =47r |Vb, Q | n a , with a representing the L or R leads, and the above equation can be rewritten as 



E< n ,(t 1 ,i a ) = ir a ff(t 1 )fl(t 2 ) 



— — exp[-ie a (ti - t 2 )\ 



x / da;'A n exp ( _(f_f^L | 1T( 



x / dxA n i exp 



P 

(x - x y 



En' 



X - X 
I 



X — Xq 



V e** A*exp(- — 



V e^ x A* exp 



{xf 



Which can be simplified 54] by carrying out the integrals and the final result is written as 



(A16) 



„ , ,_ ir ^^\.=r + V¥)\(i±* )Y 



ex p[^ + (2E)2] ± 



= S<(T)S<(t 1) t a ) = S<(t 1 ,t 2 ) 



— - exp[-i£ Q (*i - i 2 )], 



where 7 = -^L, x = ^f-, and the — , + signs stands for the left and the right leads respectively. 



(A17) 
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